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The two-loop self-energy correction is evaluated to all orders in Za for the ground-state Lamb shift of 
H-Iike ions with Z > 10, where Z is the nuclear charge number and a is the fine structure constant. The 
results obtained are compared with the analytical values for the Za-expansion coefficients. An extrapolation 
of the all-order numerical results to Z = 1 is presented and implications of our calculation for the hydrogen 
Lamb shift are discussed. 



Hydrogen and hydrogen-like ions are the simplest 
atomic systems, whose studies inspired the creation 
■and development of modern quantum electrodynam- 
|ics (QED). Despite their apparent simplicity, two-body 
■atomic systems continue to challenge physicists after 
.more than a century of research. On the experimental 
"side, the absolute frequency of the ls-2s transition in 
.atomic hydrogen has lately been measured to 1.8 parts 
^ 'in 10""^* PP, which represents an improvement of accu- 
f-^ Tacy by 4 orders of magnitude achieved during the last 
I ^twenty years. On the theoretical side, the largest er- 
^^Tor for the ls-2s transition energy stems presently from 
the uncertainty in the experimental value for the pro- 
'ton charge radius |2|, even with the most recent re- 
.analysis of electron-proton scattering data Apart 
Trom that, the major theoretical uncertainty comes from 
the two-loop QED effects. Important progress achieved 
'recently in investigations of two-loop corrections to or- 
der a'^{Za)^ HEI resulted in a much better theoreti- 
cal understanding of the hydrogen Lamb shift. Conse- 
quently, in the latest adjustment of fundamental con- 
stants , the value for the proton charge radius (which 
enters as a parameter into the determination of the Ryd- 
berg constant) was obtained by comparing spectroscopic 
data for the ls-2s transition in hydrogen with the cor- 
responding theoretical prediction. 

The results of refs. 01 E] were obtained within the 
traditional approach based on a semi-analytic expansion 
of binding QED corrections in Za and ln[(Za)^^] [Z is 
the nuclear charge number and a is the fine structure 
constant). A peculiar feature of the two- loop effects 
(first of all, the two-loop self-energy correction) is that 
this expansion converges very slowly even in case of hy- 
drogen. Numerical values of the expansion coefficients 
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are large and tend to grow with increase of the order, 
which makes estimations of uncertainties due to higher- 
order effects rather problematic. Moreover, complexity 
of calculations of the Za-expansion coefficients grows 
drastically with the increase of their order, so that the 
derivation of the complete contribution to the next or- 
der a'^{Zay does not seem feasible in the near future. 

An alternative approach is to perform the investi- 
gation non-perturbatively in the parameter Za. For 
one-loop corrections, such calculations extended over 
three decades [T] |H1 El ^| ■ A non-perturbative evalu- 
ation of two-loop QED effects is a much more demand- 
ing problem. The dominant contribution (especially in 
the low-Z region) stems from the two-loop self-energy 
correction, which is considered to be the most problem- 
atic two-loop QED effect. A calculation of this correc- 
tion to all orders in Za was started in refs. ^2 El and 
completed in a series of our investigations ^1 114L [T51 . 
Up to now, the calculation of the two-loop self-energy 
correction was carried out for ions with Z > 40 only. 
Large numerical cancellations growing fast when Z was 
decreased prevented us from calculating this correction 
for lower- Z ions and from drawing any definite conclu- 
sions about agreement of our results with the known 
Za-expansion terms. The goal of the present investi- 
gation is to perform a calculation of the two-loop self- 
energy correction for ions with Z > 10 and to compare 
our non-perturbative treatment with the previous inves- 
tigations within the Za-expansion approach. 

We start with summarizing the results available for 
the Za expansion of this correction for the ground state 
of H-like ions. The corresponding energy shift is conve- 
niently represented in terms of a dimensionless function 
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F{Za) (relativistic units h—c=l are used throughout 
the paper), 



AE = m{-) ^^F{Za). 



(1) 



where n is the principal quantum number and the Za 
expansion of the function F reads 



F{Za) = S40 + {Za)B:M + {Zaf 
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^L^B62+LBei + Gi,.o.{Za) 



(2) 



where L — ln[(Za)^^] and Gh.o. is the higher-order 
remainder whose expansion starts with a constant, 
Gh.o. {Za) = Beo + Za (...)• The results presently avail- 
able for the expansion coefficients (see refs. and 
references therein) are: 



Bio{ns) = 1.409 244... , 
B^oins) = -24.2668(31) , 
Besins) = -8/27, 
S62(ls) = 16/27- (16/9) ln2, 
Beiils) = 49.838317, 
Beo{ls) = -61.6(9.2). 
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Contrary to the calculations summarized above, our 
present consideration does not rely on the Za expan- 
sion. The working frame is the Furry picture, where the 
interaction of the electron with the nucleus is taken into 
account to all orders right from the beginning. The price 
to pay is that we have to deal with the bound-election 
propagators, whose structure is much more elaborate 
than that of the /ree-electron propagator. The detailed 
analysis of the two-loop self-energy correction was pre- 
sented in our previous investigation JHl and will not be 
repeated here. In this Letter we just mention the major 
problems to be tackled in order to obtain reliable nu- 
merical results in the \ow-Z region and sketch the ways 
of their solutions. (To note, numerical cancellations for 
the lowest value of Z considered here, Z = 10, amount 
to 4 orders of magnitude.) 

One of the main factors that influences the numerical 
accuracy of the final result is its dependence on the size 
of the basis set used in the evaluation of the P term (see 
ref. for details). In the present work, a new dual- 
kinetic-balance basis set was employed, which con- 
siderably improved the convergence in the I0W-.Z region. 
In addition, the actual calculation was repeated several 
times with increasing numbers of basis functions N and 
then an extrapolation N ^ 00 was performed. We note, 
however, that despite the improvement achieved, the de- 
pendence of the P term on the basis set is still one of 



the major sources of the uncertainty. Another difficulty 
in the evaluation of the P term was to properly control 
the accuracy of numerical integrations over momentum 
variables. It is associated with a significant contribution 
coming from the region of very large momenta, where 
the numerical Green function is not smooth enough, due 
to restrictions of a finite-basis-set representation. The 
problem has been handled by introducing a set of sub- 
tractions that have the same behavior for large momenta 
as the original integrand but are easier to evaluate nu- 
merically, and by employing very fine grids for numer- 
ical momentum integrations. The calculation of the M 
term (see ref. |15| for details) was carried out employing 
the contour Clh for the integrations over the virtual- 
photon energies, which is much more suitable for the 
numerical evaluation in the low-Z region than the inte- 
gration simply along the imaginary axis. In addition, a 
contribution containing the dominant part of the spu- 
rious behavior in the low-Z region was separated from 
the M term and calculated separately. It involves only 
one infinite partial-wave expansion, which makes its nu- 
merical computation easier as compared to the full M 
term. 

In the table, we present the numerical results ob- 
tained for the two-loop self-energy correction for the 
ground state of H-like ions with Z > 10. As compared to 
our previous investigations |14II15| . new calculations for 
Z = 10-30 were performed and the numerical accuracy 
for Z = 40, 50, and 60 was improved. The numerical 
values for .Z > 70 are taken from our previous investiga- 
tion jl5| . In fig.n our non-perturbative (in Za) results 
are compared with the contribution of the known coeffi- 
cients of the Za expansion, as given by eqs. (PJ-®- We 
observe that the numerical results behave smoothly and 
tend to approach the known analytical value at Z — 0. 

In order to perform a more detailed comparison with 
the Za-expansion calculations, we separate the higher- 
order remainder G'^ °' (Za) [defined by eq. 101] from our 
non-perturbative results, with the corresponding plot 
presented in fig. |21 As can be seen from the figure, a 
naively estimated limit of our numerical values at Z = 
is about twice as large as the analytical result (jSJ for 
the coefficient Bqq. However, there is an indication ^21 
that the analytical result of ref. |4j for the coefficient 
Bqi is incomplete. In this case, the numerical data for 
Qii.o. p\Q^^Q(^ fig. |21contain an admixture of the loga- 
rithmic contribution to the leading order and, therefore, 
are not bound at Z = 0. Until analytical calculations of 
Bqi have been finished, we cannot draw any conclusions 
concerning agreement with the existing result for Bqq. 

It is possible, however, to extrapolate our numerical 
results for G^^ ° {Za) to Z ~ \ and to obtain an approx- 
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imate value for the corresponding higher-order contri- 
bution for the hydrogen Lamb shift. For such extrapo- 
lation, we use the procedure first employed in ref. |18| 
and recently described in detail in ref. (19) . The approx- 
imate value for G^ ° {la) is obtained in two steps. First 
we apply an (exact) linear fit to each pair of two consec- 
utive points from our data set and store the resulting 
value at Z = 1 as a function of the average abscissa 
of the two points of the original set. (The points with 
Z = 10 and 15 are not employed for the extrapolation 
because of their large error bars.) Then, we perform a 
global linear or quadratic least-squares fit to the set of 
data obtained and take the fitted value at Z = 1 as a 
final result. 

In order to test this extrapolation scheme and to 
check the consistency of our numerical data with the 
first Za-expansion coefficients, we consider the function 

FiZa) = ^^^^^p^ = B,o + {Za)i...). (9) 
Za 

The extrapolation procedure described above yields a 
result that reproduces the analytical value for the B^q 
coefficient within the 15% accuracy. Applying the same 
extrapolation scheme to the higher-order remainder, we 
obtain 

Gh°(la) = -127±30%. (10) 

The error bars indicated are obtained by applying the 
extrapolation procedure to the function G^^°■{Za) -\- 
10 In(Za)^^, which is supposed to account for the prob- 
able incompleteness of the present result for the Bqi 
coefficient. 

The result (|10|l significantly alters the previous pre- 
diction for the higher-order remainder [S], which reads 
(in our notations) 

G''°(la;old) = -61.6±15%. (11) 

The difference between ()10|l and (|llll leads to a shift of 
the latest prediction for the hydrogen Is Lamb shift 5 
by 7 kHz, with the result 

i/(Ls) = 8172 804(32)(4)kHz, (12) 

where the first error stems from the current uncer- 
tainty of the proton charge radius and the second one 
is the theoretical uncertainty corresponding to the one 
of eq. I|10(l . We note that significant progress in the de- 
termination of the proton radius is anticipated in the 
near future from the experiment on the muonic hydro- 
gen, which is currently being pursued in Paul Scherrer 
Institute PU]. 



To sum up, we have performed a non-perturbative 
(in Za) calculation of the two-loop self-energy correc- 
tion that extends our previous evaluation to the region 
Z > 10. The numerical results obtained agree well with 
the first two terms of the Za expansion. A certain dis- 
agreement is found with the analytical results to order 
a^(Za)^, which could possibly be associated with in- 
completeness of the present value for the Bqi coefficient. 
An extrapolation of the numerical data to Z — \ yields 
a result that alters the theoretical value for the the hy- 
drogen Is Lamb shift by 7 kHz. 
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The two-loop self-energy correction evaluated to all or- 
ders in Za (dots) together with the contribution of 
all known terms of the Za expansion (dashed line), 
in terms of F{Za). The cross indicates the analytical 
value of this correction at Z = 0. 
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Higher-order remainder G'^'°'{Za) defined by eq. 
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Individual contributions to the two-loop self-energy correction expressed in terms of F(Za). 



z 


LAL 


F term 


P term 


M term 


Total 


10 


-0.3577 


822.14(2) 


-721.34(12) 


-100.19(10) 


0.25(16) 


15 


-0.4951 


292.902(13) 


-235.205(70) 


-57.366(48) 


-0.164(85) 


20 


-0.6015 


136.911(7) 


-102.026(55) 


-34.764(16) 


-0.481(58) 


30 


-0.7565 


44.729(3) 


-29.410(25) 


-15.465(5) 


-0.903(26) 


40 


-0.8711 


19.505(3) 


-11.575(30) 


-8.253(5) 


-1.194(31) 








-11.41(15)'^ 


-8.27(18)" 


-1.05(23)" 


50 


-0.9734 


10.025(2) 


-5.488(26) 


-5.001(3) 


-1.437(26) 








-5.41(8)" 


-4.99(6)" 


-1.34(10)" 


60 


-1.082 


5.723(1) 


-2.970(18) 


-3.341(2) 


-1.670(18) 








-2.93(4)'^ 


-3.342(21)" 


-1.63(4)" 


70 


-1.216 


3.497(1) 


-1.757(25) 


-2.412(11) 


-1.888(27) 


83 


-1.466 


1.938 


-1.057(13) 


-1.764(4) 


-2.349(14) 


92 


-1.734 


1.276 


-0.812(10) 


-1.513(3) 


-2.783(10) 


100 


-2.099 


0.825 


-0.723(7) 


-1.384(3) 


-3.381(8) 


" ref. [ig. 



